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Abstract 

We phenomenologically investigate stick-slip motion of sheared granular layers. Our phenomenol- 
ogy is constructed in the context of nucleation-and-growth of the fluidized area which is triggered 
by collapsing of stress chains. Based on this picture, we give a simple friction model by introducing 
the degree of the fluidization. It is found that the present model can successfully reproduce major 
features of the experimental results reported by Nasuno et al. (Phys. Rev. E 58, 2161 (1998)) 
with quantitatively good agreement. 

PACS numbers: 45.70.-n, 45.70.Vn 
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Granular friction is one of the central topics of granular science However, our 

understanding is still not enough both in microscopic level (friction among the individual 
particles) and in macroscopic level. 

Previously, Nasuno and his co-workers reported the detailed experimental studies of stick- 
slip motion observed in sheared granular layers [2j . Their experimental method is very simple 
and is similar to those used in studies of solid-on-solid friction ^, ^ 0]. That is, the cover 
plate with the mass M is placed on the granular layer, and then it is pushed by a spring 
(spring constant k) which is driven at a constant speed V. Despite this simplicity, the 
obtained results reveal many remarkable things which are distinguishable from earlier works 
on friction phenomena in several points. Here, lets us summarize main results of their 
experiment: (i) By the fast and sensitive measurements of the cover plate on the granular 
layers, they could determine the friction force as a function of the instantaneous sliding 
velocity. The resultant friction force is found to be a multivalued function of the sliding 
velocity during the slip event. Then characteristic hysteresis is observed, (ii) As the driving 
speed V is increased, the stick-slip motion gradually changes into the oscillatory motion, 
where the inertia of the cover plate becomes dominant. Furthermore, at a critical driving 
speed Vc the transition to the steady sliding motion occurs. For small k the transition 
exhibits pronounced discontinuity, and large fluctuations can be observed in the vicinity of 
the transition point, (iii) In the stick-slip regime, the observed instantaneous velocity during 
the slip event exhibits almost universal behavior for different driving speed V when k and 
M are fixed, (iv) Their measurements of vertical motion of the cover plate indicates that 
the dilatancy might play a crucial role in the granular friction. Above all, even within the 
sticking interval, a small creep occurs due to the localized microscopic rearrangement of the 
granular particles. From their detailed observation, they ascribed the observed stick-slip 
motion to repeated "fiuidization" and "solidification" of the granular layers. (However, the 
term "fiuidization" might not be adequate because in the experiment the cover plate moves 
at most for a distance of a few particles.) 

n 

Motivated by the experiments by Nasuno et al. p|, several researchers have proposed 
simple models to understand the experimental results 3,0]. Surely, their works succeeded 
in reproducing some aspects of the experiment, but, as a whole, further consideration seems 
to be needed (of . In this paper we also propose a simple phenomenology for sheared granular 
layers. In order to describe an internal state of the granular layer we introduce the degree 
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of the fluidization 6, similarly as Refs. [7|, |10|, but in the different context. Our model is 



constructed by interpreting the fluidization of the granular layers as nucleation-and-growth 
of the fluidized area which is triggered by the collapsing of stress chains. As shown in the 
following, the present model can successfully reproduce major features of the experimental 
results reported by Nasuno et al. 3] with quantitatively good agreement. Because of the 
considerable difficulty from the first principle approach, there is still no reliable general 
theory which can explain extensive aspects of granular dynamics. So, we believe that the 
phenomenological approach is useful for the present stage of granular physics. 

We now introduce a simple friction model for the sheared granular layers. Here, using a 
schematic illustration shown in Fig. ^ we shall explain our phenomenology as follows, (i) 
During the sticking period, the friction force is balanced with the applied force from the 
pulled spring. In the meantime, the stress chains develop in the granular layer so that they 
sustain the increased applied force, which should be a microscopic origin of increasing of the 
friction force. Eventually, the system reaches to the state exhibiting the maximum static 
friction F^aax (Fig- E^a)), where the granular layer can no longer sustain more external force, 
(ii) When the applied force exceeds the maximum static friction -Fmax, the cover plate begins 
to slide. Simultaneously, the collapsing of the stress chains should trigger the fiuidization of 
the granular layers, giving rise to pronounced decrease of the friction force. It is known that 
the network of the stress chains spreads heterogeneously in the granular layer, resulting in the 
heterogeneous stress transmission. Instead of imaging that the fiuidization does not happen 
uniformly, hence, we now ascribe it to nucleation-and-growth of the fiuidized area (Fig. 
^b)). In this paper "fiuidized area" is interpreted as a region in which the sustained stress 
is smaller than that in the region with fully developed stress chains. This is due to the break- 
down of the stress chains. Since the fiuidization is a local process at the fiuidization front, 
it is assumed that the fiuidization fronts propagate with a constant velocity independently 
of the cover plate, (iii) In a short time, the whole area becomes fiuidized (Fig. E^c)), and 
then the friction force exhibits almost saturated value, (iv) Finally, the cover plate ceases 
to move, and then the stick regime recovers again. In the experiment j^, this re-sticking 
process suddenly occurs with pronounced decrease of the friction force. 

We are at the position to construct our friction model based on the above mentioned 
picture. To begin with, we introduce the degree of the fiuidization 9 with the aid of the simple 
picture of phase change, instead of considering the microscopic dynamics of the granular 
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particles. Let us imagine, for an example, a magnetic system in an external magnetic field, 
where constituent spins point to the preferable direction. If the external field is turned 
over, the whole region of the system becomes unstable, and then stable phase, with spins 
pointed to the opposite direction, nucleates and grows. This nucleation-and-growth for non- 
conserved order parameter system has been known as phase change (or phase switching) and 
is well described by Avrami-Kolmogolov theory jlll . Il2[ llS^ . The velocity of propagating 
front of new phase is a constant, so that in the early stage of growth process, the volume 
fraction of the new phase grows as f^^^, with d being a spatial dimension, in the case of 
constant nucleation rate. On the other hand, the volume fraction grows as in the case 
where there are only latent nuclei at the beginning but no new nucleation during the growth 
regime. In this paper, using this analogous situation, the fiuidization is interpreted as a 
2-dimensional phase change with latent nuclei at the beginning. Thus, 6{t) is introduced as 
the volume fraction of the fiuidized area, and is assumed to be given by the following simple 
form 

9{t) = l-exp{-A{t)), (1) 

where A{t) is the so-called Avrami's extended volume fraction 12, 13|, which is given 
by A{t) = (ar)^ in the present case, where is some characteristic time relating to the 
microscopic mechanism of the fiuidization. The saturated state is realized (Fig. ^^c)) when 
ar ^ 1. Here r = t — nT, where n is a number, and T is the period of stick-slip event, 
respectively. We set t = just at the beginning of the first slip event. The exponent 2 
may be replaced by some fractal number reflecting the complicated structure of the growing 
front. 

An equation of motion of the cover plate is given by 

Mx = kSx - F, (2) 

where 6x = Vt — x is the displacement from the natural length of the spring. When the 
fiuidization proceeds as in Fig. ^b), the friction becomes smaller for the part of the fiuidized 
area. Therefore, we propose the following dynamic friction force in terms of 6 

Fd = Fo + Fi[l-e{t)] + bx, (3) 

where Fi[l — 6{t)] represents the decreasing part of the friction owing to collapsing of the 
stress chains, and Fq + Fi = Fmax- Practically, in an usual solid-on-solid friction, it is well 



known that the friction force is proportional to the real contact area rather than to the 
apparent contact area [l^. Eq. (jS)) can be understood as its modified version. Here we 
add the velocity dependent friction force bx, which is needed to describe correct transition 
behavior. It must be mentioned that the present model as well as the models in Refs. 0, 
cannot describe the re-sticking process, and therefore Eq. (jH)) is valid only for x > 0. So, in 
the sticking period we must use the static friction as Fg = kSx. 

We make Eq. Q a dimensionless form by using the variables, t = ^JkjMt^ x = {k/Mg)x, 
and F = F/Mg. Then, the equation of motion (0) becomes 

X = at — X — F, (4) 

where a = V/g^Jk/M, and the friction force is written as 

F, = Fo + Fi[l-e(t)] + ^l (£>0), 

Fs = at -X {at- X < Fmax)- 

In the following numerical analysis, we set h = 0.392(kg/s), Fq = 0.44, and Fi = 0.22, which 



are estimated from the experimental results, a = a is an adjustable parameter and 

we set a = 146.99 (1/s). 

Eq. (jU with Eq. (jSJ can be numerically integrated in a light manner. The obtained 
typical stick-slip motion is shown in Fig. |21 where the numerically calculated (a) deflection 
6x(t), (b) position x(t), and (c) instantaneous velocity x{t) are presented, respectively. In 
the following, we shall consider our numerical results in some detail. Fig. |H1 represents in- 
stantaneous velocitv for various V at k = 134.7(N/ni) and M = 10.90(g), which corresponds 
to Fig. 11 in Ref. [3]. Similarly to the experiment [j], we find that the instantaneous veloc- 
ity is hardly changed against V. This is entirely due to our choice of the form of F. That 
is, our friction model is intrinsically insensitive to the driving speed V. Now, it must be 
recalled that we have constructed our phenomenology imaging that the fiuidization occurs 
as nucleation-and-growth. Then, the degree of the fiuidization 6 has been introduced as 
Eq. (^. There, we have implicitly assumed that the motion of the cover plate is important 
only as a trigger of the break-down of the stress chain, and then the resultant collapsing 
dynamics of the stress chains proceeds almost independently of the motion of the cover 
plate. The experimental result might support our assumption indirectly. From Eq. 
we can calculate the friction force as a function of x. The result is shown in FigHJ where 
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the filled circles represent experimental results P]. In the early stage of the slipping event 
{Vk/ Fio^ <^t <^\/a), the velocity x is well approximated by 

3M ^ ^ 

Hence, the friction force decreases as 

F = Fmax + &i-(3M/F\«)ii;i. (7) 

Though the above expressions depend on the detailed structure of the fluidization, Eq.(IT]), 
these behavior in the early stage of the slip event is also well consistent to the experimental 
results P|. 

Next, we present numerical results for the dynamic transition behavior observed in the 
experiment [2|. As increasing the driving speed V, the inertia of the cover plate becomes 
crucial. As a result, the stick-slip motion becomes smoother, and then the oscillatory motion 
is observed (Fig. Efb)). Moreover, above a certain critical value Vc, the oscillatory motion 
changes into the steady sliding motion (Fig. El^c)). The period and the amplitude of stick- 
slip motion against the driving velocity V are shown in FiglHl in which the dotted line ceases 
at the point where the stick-slip behavior vanishes. In an usual stick-slip regime, the period 
behaves as and the amplitude of 6x{t) is almost constant. But, as V increases, the 
modulation from such behavior becomes pronounced due to the dominance of the inertia of 
the cover plate. Comparing with the experimental result by Nasuno et ai, we find that our 
obtained transition behavior shows good agreement with quantitative sufficiency. 

In this paper we have phenomenologically investigated stick-slip friction of sheared granu- 
lar layers reported by Nasuno et al. [3]- Our model constructed in the context of nucleation- 
and-growth of the fluidization can reproduce major aspects of the experimental results. We 
consider that the present idea might be also useful to the investigation of the dynamics of 
earthquake. However, our model is still insufficient to the complete understanding of the 
experimental works, (i) We can say nothing about the re-sticking process. In the experiment 
P], the remarkable decrease in the friction force is observed just before re-sticking. At this 
time, the dilation of the granular layers ends. This observation indicates that the vertical 
motion of the granular layers might play a crucial role to the re-sticking, (ii) Our present 
approach is mean-field like, so that the fluctuation effect is not included. In the experiment 
jsl, such an effect surely leads to intriguing chaotic stick-slip motion for very low V and 



large k. These problems are left for future works 
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(a) (b) (c) 

FIG. 1: Schematic illustration of our phenomenology seen from the top of the granular layer, (a) At 
the maximuin friction -^^max" 

(b) When the applied force exceeds -Pmax, the granular layer begins to 
fiuidize. Our phenomenology assumes this fluidization proceeds as nucleation-and-growth, and the 
nuclei is attributed to the points at which the break-down of stress chains initially take place. The 
fluidized areas represented by dots dynamically coexist with the non-fluidized areas represented by 
diagonal lines, (c) The whole region becomes fluidized. 
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FIG. 2: The numerically calculated (a) deflection Sx{t), (b) position x{t), and (c) instantaneous 
velocity x{t), for k = 134.7(N/m), M = 10.9(g), and y=113.33(^m/s). Here the parameters used 
are the same as those in Fig. 3 of Ref . . 
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FIG. 3: Instantaneous velocity of the cover plate x for various value of y at = 135 (N/m) and 
M = 10.90(g). The parameters used are the same as those in Fig. 11 of Ref. (31] • 
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FIG. 4: The normalized friction force F/Mg as a function of instantaneous velocity x for k 
134.7(N/m), M = 10.9(g), and V = 113.33(//m/s). The filled circles represent the experimental 
results, which are taken from Fig. 13 of Ref. 
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FIG. 5: Deflection 6x{t) for various V at k = 1077(N/m) and M =10.9(g). (a) Typical stick-slip 
motion is observed at V = 5.67(;um/s). (b) Oscillatory stick-slip behavior at V = 5.67(mm/s). (c) 
Steady-sliding motion at V = 12.4(mm/s), where we do not show the transient damped oscillatory 
motion observed for t < 0.5(s). The parameters used here is the same as those in Fig. 7 of Ref. 
Pi, except for (c), where we use slightly different value of V. 
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FIG. 6: (a) Amplitude of the deflection 5x of stick-slip motion as a function of V for various k. 
(b) Period of stick-slip motion as a function of V for various k. The open marks are taken from 
Fig. 8 in Ref. Q], where the data for k = 135(N/m) is not available for 6x. 
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